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Abstract—In order to satisfy a relatively high resolution for the retrieval of snow water equivalent, an X/Ku-band dual-frequency 
full-polarized SCATterometer (DFPSCAT) on-board Water Cycle Observation Mission (WCOM) satellite is designed for 
high-resolution observations. However, given the situation that i) the conically rotation rate of antenna is relatively fast (19 r/min), ii) the 
swath width is larger than 1000 km, iii) day or night observation capabilities are required, the method called “Rotating Azimuth 
Doppler Discrimination” is proposed, which can satisfy the resolution requirement and real-time processing. Considering the 
complexity of the system's design and the improvement of azimuth resolution capability, a burst pulsing scheme is addressed to satisfy 
the numbers of azimuth sampling. The simulation model is used to analyze the feasibility of azimuth discrimination method based on 
geometry and system parameters. It is shown that the achievable azimuth resolution is about 2-5km far end of the swath and only 5km 
near end of the swath. The results show that Kpc is <0.4 when snow thickness exceeds 5 m and about 0.75 when light snow. 


Index Terms—Scatterometer, Rotating Antenna, Doppler Discrimination, burst pulse scheme 


I. INTRODUCTION 


WCOM proposes a DFPSCAT using the dual frequency (X-band 9.6 GHz and Ku-band 17 GHz) and multi-polarization for 
mapping of snow water equivalent and freeze-thaw processes [1, 2], aiming to provide higher accuracy and consistent 


measurements of key elements of water cycle from space outside of the Earth's atmosphere, including soil moisture, ocean salinity, 
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freeze-thaw, snow water equivalent, etc. Space-borne scatterometer has become of great benefit to the numerical weather 
predictions and climate monitoring in the past 40 years, and the design and development of a scatterometer for different 
applications has been applied. The measurement accuracy and performance of snow water equivalent is related to the resolution. In 
order to take into account the local vegetable variability, the resolution should be degraded to a higher resolution of 1km [3]. 
However, such a high resolution is hard to achieve for current scatterometer instrument, which has an inherent spatial resolution on 
the order of 10-50 km [4, 5]. So, considering the implementation capability of the DFPSCAT on-board WCOM, a spatial resolution 
of 2-5 km is proposed [2]. 

Unlike synthetic aperture radar (SAR), real aperture radar (RAR) provides measurements in real-time processing. However, 
especially when the antenna is scanning fast, the accumulating time of the RAR is not long enough, so the resolution will be 
relatively low. In contrast to a fan-beam system, the pencil-beam design employs a single antenna which is conically scanning 
about the nadir axis in order to provide multiple azimuth measurements. Such a pencil-beam antenna is also easier to be installed on 
the spacecraft [6]. Based on the platform motion and the antenna scanning style of the on-board scatterometer, this system will 
achieve an improved resolution along the azimuth direction by using the Doppler discrimination method. This method is applied 
when the azimuth angle is in the range between 10 to 170 degrees [7, 8]. When the azimuth is between the intervals of 0+/-10 and 


180+/-10 degrees (the absolute value of Doppler shift is the largest and the Doppler bandwidth is near zero), the method called 


super-resolution reconstruction is applied [9]. However, that is not the focus in this paper. Therefore, a tradeoff between azimuth 
Section |] introduces the system parameters and a burst pulse scheme. In section [I], we analyze the azimuth signals, and 


introduce a rotating azimuth Doppler discrimination method. Section IV gives some simulations of point targets. Finally, the 


summary and conclusion remarks are given in Section V. 


II. INSTRUMENT DESIGN OF THE SCANNING PENCIL-BEAM SCATTEROMETER 


A. System parameters 


The DFPSCAT on-board WCOM will observe the Earth for 3 years from a sun-synchronous, near-circular orbit with an 
equator-crossing altitude near 600 km at two frequencies with four polarization (VV, HH, VH, HV) [10]. The space-borne 
scanning pencil-beam scatterometer can provide high signal-to-noise ratio (SNR) and the advantage of covering a larger area than 
SAR. The DFPSCAT instrument is a dual-frequency polarized radar scatterometer, which consists of a mechanically scanned 
radiometer with a 1.5m diameter reflector antenna, which can satisfy the requirement of high Signal Noise Ratio (SNR) and high 


gain. Some key specifications are shown in TABLE I. 
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The DFPSCAT instrument is based on a single reflector measuring up to 1.5 X 1.5 m, illuminated by two adjacently placed feeds 


that generate the multiple beams. A single rotating antenna allows a larger swath of measurements at a constant incidence angle. In 


order to get a ground swath of about 1000km, the antenna incidence angle is at least 39° , which is calculated by using the 


observation geometry, see Fig. 1. Considering the ‘along-track continuity constraint’, the rotation rate should be given by 


v 
Q > ground ( 1 ) 
X 


elevation 


where Q) is the antenna rotation rate, V 


ground 18 the ground velocity of the satellite and x, is the elevation width of the 


levation 
footprint (about 23 km in Ku-band and 28 km in X-band). The rotation rate should be at least 19 revolutions per minute (RPM) 


while considering achieving complete surface coverage over the swath. The antenna scanning loss is defined by 


n _ fe (R)8,(R)aA, 
‘scan 2 R)dA 
fe? (R) aA, > 


where g (R) and g, (R) are the beam patterns on the surface at the time transmit and receive. The integral range is the whole 


footprint. The scanning loss of a reflector antenna is represented as a function of the beam patterns on the surface at the time of 
transmit and receive, see Fig. 2. The speed of the footprint due to the satellite motion in the azimuth direction is far smaller than the 
speed caused by the scanning motion, so the scanning loss is an important consideration. Given the rapid scanning motion of the 
antenna, an analysis based on the radar ambiguity function is found to be particular useful for the pencil-beam scatterometer. Here, 
the scanning loss is a function of azimuth displacement normalized by the two-way azimuth beam-width. Because of the same 
transceiver antenna, the scanning loss should be carefully calculated. When the antenna rotates at 19 RPM, the antenna scanning 


loss is about -1.3dB. 
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Fig. 1. Geometry of the DFPSCAT instrument conically scanning 5 10 15 a 25 a 


antenna beam mapping out a swath of 1000 km at Earth’s surface. The . the rotation rate / RPM A 
antenna footprint is represented by a small oval. Fig. 2. The scanning loss vs. antenna rotation rate. 


TABLE I 


SYSTEM PARAMETERS 
X band Ku band 
Frequency 9.6 GHz 17 GHz 
Wavelength 3.10 cm 1.76 cm 
Bandwidth 2 MHz 
Orbit height 600 km 
Satellite velocity 7500 m/s 
Antenna incidence 39° 


Polarization 
Antenna 


HH, VV, HV/VH 
scanning pencil-beam 


Rotation rate 19 RPM 
Antenna dimension 1.5 mx1.5 m 
Swath width >1000 km 
Scanning loss -1.3 dB 
System loss -5 dB 
System temperature 300 K 
Peak power 200 W 150 W 
Antenna gain 42 dBi 47 dBi 
Footprint 28.0 km(ele)x 23.0 km(ele)x 
20.3 km(azi) 13.9 km(azi) 
Resolution 
(after processing) 2~5 km 


B. Radar pulse repetition frequency 


The resolution of azimuth in traditional scatterometer is equal to the azimuth footprint, which is very coarse. In order to 
unambiguously resolve the scene in azimuth, the PRF should be at least larger than the total Doppler bandwidth across the azimuth 
footprint. For Ku-band, the azimuth footprint corresponds to a Doppler bandwidth of approximately 10 kHz, and thus a minimum 
10 kHz PRF is required. Meanwhile, in order to avoid range ambiguities within the footprint, a maximum PRF of 9.46 kHz is 
required, corresponding to the beam fill time of about 0.106 ms. The resolution is only improved in one dimension for the 


traditional pulse design, so a tradeoff between range ambiguities and azimuth ambiguities should be considered [11]. One way to 


resolve this problem is to use a burst pulsing scheme [6]. In Fig. 3, the echo returns are depicted as semicircles, and Tr is the time 


interval of pulses in a transmit burst, G is the length of the transmit burst, Toe is the length of the receive burst. A multiple-pulse 


burst of length is repeated at the burst repetition interval (BRI). 


The DFPSCAT radar emits H-polarized and V-polarized pulse pairs at the burst repetition interval of approximately 250 Hz, 


which T, is 1.2 ms, and T. 


int 18 75 microseconds, which is modulated to a linear modulation frequency (LFM) chirp of 2 MHz 
bandwidth. In a transmit burst, the radar emits 16 H-polarized and 16 V-polarized pulses. For the transmit burst, pulses will adopt 
different carrier frequency of the adjacent pulse. In the process of echo returns received, pulses will be distinguished by de-chirp 
method, see Fig. 4. In a receive window, there will be several pulses with different carrier frequencies. In the process of simulation, 


the voltage of H polarization is double than V polarization. Although these pulses are overlapping in the time domain, the pulses 


can be separated by using a receive window corresponding to different carrier frequency [12]. In Fig. 4(a), there are five pulses in 
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a receive window. The V 


echo2 1S What we want in this receive window. In the Fig. 4(b), the pulses are processed by using a simple 


lowpass filter, and such processing cannot compensate the phase variation or frequency shift [13] 
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Fig. 3. Time series of radar transmit, reception, polarization. 
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Fig. 4. (a)pulses in a receive window, without filtering (b)after filtering. 


III. THE AZIMUTH RESOLUTION 
The antenna of scatterometer is not only flying at the speed of about 7550 m/s, but also rotates about the nadir axis at a stable rate. 
The roundtrip flight time for the beam pulse is approximately 5.26 ms, which can cause a relatively large rotation angle of the 
antenna (about 0.6° ). The effective area illuminated by the antenna is defined by 3 dB bandwidth. If the dimension of an antenna 
is designed large enough, the azimuth resolution can achieve several km, but a large antenna will brings great difficulties to 
implement on-board. In order to achieve global seamless coverage, in this paper a pencil-beam antenna is employed at the rate of 


19 RPM based on the along-track continuity constraint, as mentioned earlier. 


A. The analysis of azimuth resolution 


Based on the design of a burst pulsing scheme, the PRF can reach about ten kHz. As previously discussed, PRF should be higher 
than the Doppler bandwidth of -3 dB footprint. Errors caused by the satellite motion and orientation will lead to the errors for the 
Doppler centroid real-time estimation. However, this system has a high precision of the orbit information. The performance 
requirement of the attitude control error reaches 0.05 degrees, and in actually the attitude control error can reach 0.01 degrees. 
When the antenna is perpendicular to the flight direction (azimuth=90 or 270 degrees), the system gets the maximum Doppler 
bandwidth, and the iso-delay and iso-Doppler contours are locally perpendicular to the elevation and azimuth axes respectively. 
When the antenna scans to the forward or aft looking directions (azimuth=0 or 180 degrees), the iso-Doppler lines are parallel to the 
azimuth axis. In this case, the Doppler bandwidth of the azimuth footprint is almost zero. Thus it is difficult to achieve high azimuth 
resolution when the antenna scans to the forward or aft looking directions. From Fig. 5(a), the cross-track distance (CTD) is defined 


as the vertical distance between a given measurement and the nadir track. 


As the antenna is rotating, the achievable azimuth resolution Ô, is defined as: 


= Ry: A 
aT 


sat 


f (8) (83) 


where f (0...) represents the degradation in the azimuth resolution due to squint elongation effects, 0, is the azimuth angular, 


VA 


A is the wave length, R) is the slant range, V, 


sat 


is the speed of the satellite, and Ta is the effective dwell time, which is 


approximately one half of T, . The dwell time is the reciprocal proportion to the Doppler resolution. The burst pulsing scheme 


comes at the price of a factor of two or three decrease in achievable azimuth resolution, compared with traditional pulse design. 
From Fig. 5(b), when the CTD is in the range 200 to 500 km, the azimuth resolution can reach 2 km and when the CTD is about 80 
km, the azimuth resolution can only reach 5 km. When the CTD is less than 80 km, some other methods are used, which are not the 


focus in this paper [9]. 
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Fig. 5. (a) Range-Doppler geometry of rotating radar beam. (b) the azimuth resolution with CTD 


B. A rotating azimuth Doppler discrimination 

The conically scanned pencil-beam antenna will rotate as it moves along the flight direction. Meantime, the system should meet 
the requirements for real-time processing and long term data storage. So the processing algorithm must be simple and efficient, 
which suits the needs of real use on the satellite. The processing algorithm which is called “Rotating Azimuth Doppler 
Discrimination”, can meet the design requirements. 

From geometry of the observation, the echoes return from a footprint can be written as: 


_ : T—2R(tsr,x)/c : t 
s, (t,t) = ff o(r,x) reaf ) rect | 


footprint P b 


p 2 
p| jZ reira -exod jay| rt premna f) 


Cc 


(4) 


where rect( t / I ) is a rectangular window, O (r 5X ) is the backscattering coefficient of target (r PD ) 3 Jo is the transmit carrier 
frequency, T, is the width of the transmit pulse, 7 is the frequency modulation rate, C is the speed of light, 7 and t is the time 
of range and azimuth, which is called ‘fast’ time and ‘slow’ time, R (t ST x) is distance between the target and the transmitter. 


From Fig. 6 (a), the relations between the slow-time and frequency domain Doppler are described [14]. B q 18 the target Doppler 


bandwidth, T, is the dwell time of the target, k, is the Doppler rate, and k,., is the Doppler centroid varying rate, which is given 


t 


as follows: 


k _ EAO) = 2V a sin (8t) : Oo : sin(G,,,.) (5) 
rot at ~ A 


where the @,,, is the angular velocity, 0, 


is the antenna incidence angle. From Fig. 6 (a) and Fig. 6 (b) B q is the maximum at 
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target S1 and S7, and B q İs the minimum at target S4. As the antenna is scanned forward or backward, the Doppler bandwidth and 
Doppler centroid will change nonlinearly. 
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(b) 
Fig. 6. (a) The raw-data support in the slow time/frequency (b) range migration 
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Fig. 7. The block diagram of the model. 


The block diagram of the model is shown in Fig .7 [15]. The first step is to remove the carrier frequency from raw data, then the 


next step is to adopt the de-chirp method by introducing azimuth product between the signal and the following chirp signal: 
. P 2 
H, =exp[ j27 for, ]-exp| -jar (1-14) | (6) 
where fo is the carrier frequency, R; is the reference slant range, and T} =2R, / C. 


Then next three steps are to remove the residual video phase by using the range FFT, H, and range inverse FFT, where H, is 


given as follows: 


H, = exp| -jz fe l (7) 
y 


where f, is the frequency sampling in range direction. 


IV. SIMULATION EXPERIMENT 


A. Point targets response analysis 


The simulation parameters are described in TABLE I. A group of transmit burst is used, and this produces a PRF of about 13.3 
kHz. The model is used to analyze point targets firstly for the side-looking case, and then extend these results to the case of an 
arbitrary azimuth angle. The spacecraft nadir on the surface is defined as the coordinate origin. The flight direction can be defined 
as X axis and the cross track direction is Y axis. When the azimuth angle is 90° , two targets are 5 km apart, which their 
coordinates can be given as: (-2.5e3, 500e3, 0), (2.5e3, 500e3) and two targets are 2 km apart, which their coordinates can be given 


as: (-le3, 500e3), (1e3, 500e3), given in units of meter. 
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Fig 8. (a) azimuth is 90° , 5km (b) azimuth is 90° , 2km 
(c) azimuth is 60° , 5km (d) azimuth is 60° , 2km(e) azimuth is 30° , 5km (f) azimuth is 30° , 2km 
From Fig 8, the targets which are 2 km apart can be distinguished by using the model when the azimuth is in the range 60 to 90 
degrees, which is consistent with previous analysis. However, when the azimuth is 30° , the targets which are 2 km apart can not 


be separated because the Doppler resolution is slightly larger than the Doppler bandwidth between them. At the same scope of the 


azimuth, the larger the distance between the targets, the better the separation capability. The transmit burst time is 1.2 ms, 


corresponding to a scanning distance of 1.2 km. Thus the total coverage distance along the rotating direction in one burst time is 
about 15 km. When the azimuth is 30° , if the azimuth resolution is 5 km, then the scanning distance should be at least 10 km. 
Because the azimuth resolution quickly degrades near the nadir track because of the elongation effects, so when the azimuth is in 


the range 10 to 30° , the azimuth resolution should be set to about 5 km, Therefore, the azimuth resolution can reach 2 km at far 


end of swath and 5 km at near end of swath. 


B. Communication error Kpc 


The calculation of Kpc for a pencil-beam rotating scatterometer may be quite complex, especially for the combined slices [16]. 


The detected energy from several slices is summed, and a composite calibration factor is applied [17]. The normalized standard 


deviation of the Oy estimate can be expressed as: 


12 
1 2. Ly 
K p Ir 1+ —— + a (8) 
N,N. SNR SNR 
where SNR is the signal-to noise ratio, N „ is the number independent looks of the elevation, and N, is the equivalent number 


of azimuth looks. 
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| 9(a) 5km combined slice, Kpc with footprint along the cross track (b) Kpc with snow thickness using different combined slices 


The simulation data is based on the model of the dense media radiative transfer (DMRT) [18]. We choose the spatial resolution 
of a snow vector cell as 5 km x5 km and assume that the energy estimates are uncorrelated from slice to adjacent slice. From Fig. 9 
(a) and (b), the Kpc will reduce with the snow accumulation, because of the higher SNR. When snow thickness exceeds 5 m, the 
Kpc of different combined slices can reach 0.4. However the Kpc is relatively large at the shallow snow for the SNR is low. The 
Kpc depends on the combination of the slice number, the slice size and the SNR. From Fig. 9 (b), the Kpc of 3 km and 5 km 
combined slices is relatively small compared with 4 km combined slices. 

V. CONCLUSION 

In this paper, a ‘Rotating Azimuth Doppler Discrimination’ method is introduced and analyzed, which can acquire azimuth 
resolution of 2-5 km at the far end of the swath and only 5 km at the near end of the swath. Some various options considered for 
improving the azimuth resolution are compared. A key design of the burst pulse scheme is proposed which can reach 13.3 kHz 
(PRF). Because the scatterometer obtains a large swath of 1000 km in real-time and day or night, a rotating azimuth Doppler 


discriminationa technique in azimuth discrimination is employed to improve the resolution simply and efficiently. Although the 
azimuth resolution can reach a high level of 2-5 km at the near end of the swath, the normalized standard deviation of the Og 


estimate is relatively high, which may exceed 0.5. A series of simulations for the point targets have been presented to verify the 
feasibility and validity of this method. It shows that the scatterometer can obtain a high azimuth resolution by using a ‘Rotating 
Azimuth Doppler Discrimination’ method. In the future, the simulation and experiment using different parameters will be carried 


out. 
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